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ABSTRACT 

<n : 

We investigate the evolution over the last 6.3 Gyr of cosmic time (i.e., since z ~ 0.7) 
of the average X-ray properties of early-type galaxies within the Extended Chandra 
Deep Field-South (E-CDF-S). Our early-type galaxy sample includes 539 objects with 

(f) ' red-sequence colors and Sersic indices larger than n = 2.5, which were selected jointly 

from the COMBO-17 (Classifying Objects by Medium-Band Observations in 17 Filters) 
J> ■ and GEMS (Galaxy Evolution from Morphologies and SEDs) surveys. We utilize the 

deep Chandra observations over the E-CDF-S and X-ray stacking analyses to constrain 
primarily the average X-ray emission from "normal" early- type galaxies (i.e., those that 
are not dominated by luminous active galactic nuclei [AGNs]). In our analyses, we study 

\Q \ separately optically luminous (Lb ~ 10 10 ~ n Lb,®) an d faint (Lb ~ io 9 - 3 ~ 10 Lb,®) 

galaxy samples, which we expect to have soft (0.5-2.0 keV) X-ray emission dominated 

i-£h ■ by hot (~1 keV) interstellar gas and low mass X-ray binary (LMXB) populations, re- 

spectively. We detect individually 49 («9%) of our galaxies in the X-ray band, and 
classify these sources as either normal early-type galaxies (17 galaxies) or AGN candi- 
dates (32 galaxies). The AGN fraction of our optically luminous samples evolves with 
redshift in a manner consistent with the (1 + z) 3 evolution observed in other investi- 

• i-H ! gations of X-ray-selected AGNs. After removing potential AGNs from our samples, we 

find that the X-ray-to--B-band mean luminosity ratio (Lx/Lb) for optically luminous 
early-type galaxies does not evolve significantly over the redshift range z 0.0-0.7. 
This lack of X-ray evolution implies a general balance between the heating and cooling 
of the hot interstellar gas. If transient AGN activity is largely responsible for maintain- 
ing this balance, then we infer that mechanical power must be dominating the feedback 
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out to z ~ 0.7. Furthermore, in this scenario the average mechanical AGN power must 
remain roughly constant over the last half of cosmic time despite significant evolution 
in the average AGN radiative luminosity. For our optically faint early-type galaxies, we 
find suggestive evidence that Lx/Lb increases with redshift over the range z ~ 0.0-0.5; 
however, due to limited statistical constraints on both the local Lx/Lg ratio and the 
AGN contamination of our samples, we consider this result to be marginal. 

Subject headings: cosmology: observations — surveys — galaxies: normal — galaxies: 
active — galaxies: elliptical and lenticular, cD — X-rays: galaxies — X-rays: general 

1. Introduction 

Distant {z ~ 0.1-3), isolated early- type galaxies have been studied extensively at optical and 
near-IR wavelengths. The optical luminosities of early-type galaxies with rest-frame red-sequence 
colors have faded by f«l magnitude since z = 1 (e.g., Bell etal. 2004a; Mcintosh etal. 2005; Treu 
etal. 2005, 2006; van der Wei etal. 2005). Furthermore, little star-formation activity is observed 
within these galaxies out to at least z rs 0.7 (e.g., Bell et al. 2005) indicating that this observed 
evolution is consistent with the passive aging of old stars. However, the luminosity density of red- 
sequence galaxies has been observed to remain approximately constant over z ~ 0.1-0.7 (e.g., Bell 
etal. 2004a; Faber etal. 2006), and since these galaxies appear to be fading passively with cosmic 
time, this suggests that an important fraction (~l/2) of the current z = stellar-mass content of 
the early- type galaxy population has emerged since z = 1. One explanation for these observations is 
that relatively young, lower- mass galaxies continuously populate the red-sequence (e.g., van der Wei 
etal. 2004; Thomas et al. 2005) through either a downsizing mass-assembly scenario (e.g., Cowie 
etal. 1996; Bundy etal. 2005; Cimatti etal. 2006; De Lucia etal. 2006; Lee etal. 2006) and/or a 
scenario where lower-mass early-type galaxies have different merging histories than their higher- 
mass counterparts (e.g., Khochfar & Burkert 2003). At z ^ 1, distant red galaxies (DRGs; selected 
using J — K colors) and extremely red objects (EROs; R — K or I — K selected), some of which 
are likely progenitors to early-type galaxies, appear to have large star-formation rates and inferred 
dynamical masses similar to those observed for massive local ellipticals (e.g., Forster-Schreiber 
etal. 2004; McCarthy 2004; van Dokkum etal. 2004; Papovich etal. 2006). These analyses suggest 
that the majority of the stellar mass within many massive early-type galaxies formed at z ^ 1.5, 
possibly during a dust-obscured far-infrared-luminous phase (e.g., Chapman etal. 2005; Swinbank 
etal. 2006). 

1.1. X-ray Properties of Early- Type Galaxies 

X-ray investigations of local early-type galaxies have revealed that the total X-ray luminosity 
is correlated with optical luminosity (e.g., Trinchieri & Fabbiano 1985; Canizares etal. 1987; Brown 
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& Bregman 1998; O'Sullivan etal. 2001; Ellis & O'Sullivan 2006). X-ray emission from local early- 
type galaxies originates primarily from a combination of hot interstellar gas and low mass X-ray 
binaries (LMXBs), and due to their relatively low star-formation rates (SFRs) per unit stellar 
mass (typically SFR/M* « 10~ 3 -10" 1 Gyr" 1 for z « 0.0-0.7 galaxies with M* > 10 10 M ; e.g., 
Feulner etal. 2006), these galaxies are expected to have negligible contributions from high mass 
X-ray binaries (HMXBs; e.g., Grimm etal. 2002). Two fairly distinct classes of early-type galaxies, 
which can be divided by optical luminosity at Lb , C rit ~ 10 10 Lb,®, have been identified. Early- type 
galaxies with Lb ~ ^s, C rit are referred to as "optically luminous" and those with Lb ~ Lb, ait are 
referred to as "optically faint." 

Optically luminous early-type galaxies in the local universe are generally observed to have 
X-ray emission powered largely by hot (kT « 0.3-1 keV) interstellar gas at soft (0.5-2.0 keV) X-ray 
energies and LMXBs at hard (2-10 keV) X-ray energies (e.g., Matsumoto et al. 1997). Figure la il- 
lustrates the X-ray spectrum of NGC 1600 (Sivakoff et al. 2004), a representative optically luminous 
early-type galaxy; the hot-gas and LMXB components have been separated for illustrative purposes. 
The X-ray luminosities for optically luminous early-type galaxies are correlated with their optical 
luminosities following Lx <x L? B . Stellar winds generate the gas and high- velocity interactions and 
Type la supernovae heat it to X-ray temperatures (e.g., Sarazin 1997). Galaxies with larger stellar 
masses and deeper gravitational potential wells generally have higher stellar velocity dispersions 
and thus provide more significant heating of the gas. Therefore, the Lx~Lb relation is manifested 
as a result of proportionalities between the stellar mass in a particular galaxy (traced by Lb) and 
the stellar velocity dispersion of that galaxy (traced by Lx; e.g., Mahdavi & Geller 2001). 

Optically faint early-type galaxies in the local universe have soft X-ray emission characterized 
by a power-law (r « 1.8) component from LMXBs and a mild contribution from hot interstellar 
gas (i.e., a thermal plasma with kT ~ 0.1 keV). At hard X-ray energies, the emission is completely 
dominated by LMXBs. The total X-ray luminosities of optically faint early- type galaxies are 
observed to be linearly correlated with their optical luminosities (i.e., Lx oc Lb)- This linear 
correlation is thought to arise largely as a result of the number of LMXBs in a given galaxy 
being proportional to the stellar mass within that galaxy (e.g., Gilfanov 2004). Due to their low 
luminosities, few optically faint early- type galaxies have been studied in detail in the X-ray band, 
and the presently available data for these galaxies are plagued by large scatter due to small numbers 
of LMXBs per galaxy and varying contributions from hot interstellar gas (e.g., David etal. 2006). 
Therefore, much of what is known about LMXBs in these galaxies comes from studying optically 
luminous early- type galaxies that have little contamination from hot interstellar gas (i.e., X-ray faint 
early-type galaxies). Figure lb shows the X-ray spectrum of NGC 4697 (e.g., Sarazin etal. 2001), 
an X-ray faint (yet optically luminous) early-type galaxy, which is expected to be representative of 
the typical optically faint early-type galaxy X-ray spectrum. 
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1.2. Physical Properties of the X-ray Emitting Components 

The two early-type galaxy classes discussed above (optically luminous and faint) generally 
differ in their relative spectral contributions from hot interstellar gas and LMXBs. An important 
component of the X-ray emission from early-type galaxies, most notably in optically luminous 
sources, is the hot interstellar gas. Detailed observations of the gas have shown that it radiates 
strongly but does not appear to be cooling as expected (e.g., Xu etal. 2002; § 7.1 of Mathews & 
Brighenti 2003 and references therein; Tamura etal. 2003; Bregman etal. 2005). The observed gas 
temperatures and densities of the central regions of optically luminous early- type galaxies (e.g., 
Fukazawa etal. 2006) imply short radiative cooling times of ~10 8 yr on average. Simple cooling- 
flow models, which include heating from stellar winds and Type la supernovae, vastly overpredict 
the amount of cooled gas observed in the central regions of these galaxies (see, e.g., § 8 of Mathews 
& Brighenti 2003); this suggests some additional form of heating must be present. Attempts to 
mitigate this problem by including additional heating components in these models, such as transient 
heating from active galactic nuclei (AGNs), have been unsuccessful in reproducing the observed 
properties of the gas (e.g., the radial profiles of electron temperature, density, and metallicity; 
see, e.g., Brighenti & Mathews 2003). Despite these modelling difficulties, both observational and 
theoretical studies have suggested that energy feedback from transient AGNs must be important 
(e.g., Binney & Tabor 1995; Ciotti & Ostriker 1997, 2001; Birzan etal. 2004; Best etal. 2005, 2006; 
Churazov et al. 2005; Di Matteo et al. 2005; Scannapieco et al. 2005; Bower et al. 2006; Brighenti 
Sz Mathews 2006; Croton etal. 2006; Hopkins etal. 2006; McNamara etal. 2006); these studies 
maintain that AGN feedback plays a crucial role in regulating the growth of early-type galaxies and 
their central supermassive black holes. Observations of giant elliptical galaxies in the local universe 
have revealed significant X-ray cavities filled with radio lobes that are inferred to have been driven by 
nuclear AGN outbursts (e.g., Churazov et al. 2001, 2002; Jones et al. 2002, 2005; Forman et al. 2005; 
Nulsen etal. 2005a, 2005b; Allen etal. 2006); these observations suggest significant energy transfer 
and recurrence timescales much shorter than the average cooling time. One of the goals of this 
paper is to characterize the evolution of the hot X-ray-emitting gas in normal early-type galaxies 
and to place constraints on the role of feedback by transient AGNs. 

LMXBs provide another important contribution to the X-ray emission from early-type galaxies 
(both optically luminous and faint), and many investigations have been dedicated to understand- 
ing their nature. Observations in the optical band, most notably from the Hubble Space Telescope 
(HST), have shown that many of the LMXBs in early-type galaxies are coincident with globular 
clusters (e.g., Sarazin et al. 2003), and the number of LMXBs observed in a given galaxy appears to 
be linearly correlated with the globular-cluster specific frequency, Sn (i-e., the number of globular 
clusters per unit optical luminosity; e.g., White etal. 2002; Kim & Fabbiano 2004; Irwin 2005; 
Juett 2005). White etal. (2002) examined several optically luminous, LMXB-dominated early- 
type galaxies with stellar populations having ages ranging from ss5-13 Gyr (as measured via spec- 
troscopy), and found that their X-ray-to-optical luminosity ratios (Lx[LMXB]/L op t) have no de- 
pendence on galaxy age; Lx (LMXB) /L opt was observed to be linearly proportional to Sn. Again us- 
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ing optically luminous early-type galaxies, Irwin (2005) found that details of the Lx(LMXB)/L opt -SW 
relation imply that a significant number of LMXBs were formed in the fields of the early-type 
galaxies and that the field LMXB contribution to Lx(LMXB)/L opt is most significant, and even 
dominates over the globular cluster LMXB emission, for galaxies with smaller values of Sn- In 
the optically luminous galaxies studied by White etal. (2002) and Irwin (2005) Sn is large, and 
LMXBs are mainly formed in the high-density central regions of globular clusters through tidal 
interactions between normal stars and either neutron stars or black holes. In optically faint low- 
mass systems, however, Sn is typically small (e.g., Ashman & Zepf 1998), and it is therefore likely 
that LMXBs that formed within the galactic fields of these systems dominate Lx(LMXB)/L opt 
(see, e.g., Fig. 3 of Irwin 2005); these LMXBs are thought to have formed via the evolution of 
primordial binaries on timescales of ~1-10 Gyr following a star-formation event (e.g., Verbunt & 
van den Heuvel 1995). Furthermore, if the downsizing mass-assembly scenario discussed above is 
correct, field LMXBs within optically faint galaxies could evolve significantly over cosmic time, as 
LMXB populations emerge in the wake of relatively recent bursts of star formation and eventually 
fade after ~1 Gyr (e.g., Ghosh & White 2001). Therefore, another goal of this paper is to measure 
the X-ray properties of normal optically faint early-type galaxies over a significant range of redshift 
to place constraints on the evolution of the LMXB activity in these systems. 



1.3. Observations of the X-ray Emission from Distant Early- Type Galaxies 

At present, there have been few investigations of the redshift evolution of the X-ray prop- 
erties of isolated early-type galaxies, which contrasts with the case for late-type galaxies (e.g., 
Hornschemeier et al. 2002; Norman et al. 2004; Laird et al. 2005; Kim et al. 2006; Lehmer et al. 2006). 
Brand et al. (2005) used X-ray stacking analyses (see also, e.g., Brandt et al. 2001; Nandra et al. 2002; 
Lehmer et al. 2005a) to investigate the evolution of optically luminous red galaxies within the Bootes 
field of the NOAO Deep Wide-Field Survey (NDWFS) over the redshift range z ~ 0.3-0.9. The 
Chandra observations used in these analyses were «5 ks in duration (e.g., Murray etal. 2005) and 
covered an area of deg 2 . Brand etal. (2005) found that the average X-ray luminosities of 

these galaxies mildly increase with redshift, which is primarily due to a rise in AGN activity with 
redshift. Due to its relatively large sample size (3316 red galaxies), the Brand etal. (2005) inves- 
tigation provided useful statistical constraints on the X-ray activity from distant powerful AGNs; 
however, the X-ray sensitivity of this study was too low (i.e., the luminosity detection limit was 
Lx ~ 10 43 - 2 ergs s _1 at z « 0.7) to distinguish effectively galactic X-ray emission (i.e., hot in- 
terstellar gas and LMXBs) from that of luminous AGNs. Deep X-ray surveys, which utilize long 
X-ray exposures with Chandra or XMM-Newton (see, e.g., Brandt & Hasinger 2005 for a review), 
are required to investigate the X-ray emission from distant "normal" (i.e., not predominantly pow- 
ered by luminous AGNs) early-type galaxies. An ideally suited field for such an investigation is 
the Extended Chandra Deep Field-South (E-CDF-S). The E-CDF-S is a «0.3 deg 2 , multiwave- 
length survey field composed of the central «1 Ms Chandra Deep Field-South (CDF-S; Giac- 
coni et al. 2002; Alexander et al. 2003) and four flanking, contiguous «250 ks Chandra observations 
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(Lehmer et al. 2005b). The 0.5-2.0 keV sensitivity limits (for detecting an individual unresolved 
X-ray point source) over the E-CDF-S are «5.2 x 10~ 17 ergs cm~ 2 s _1 in the most sensitive regions 
and «3.0x 10~ 16 ergs cm~ 2 s _1 over the majority of the field; these levels correspond to 0.5-2.0 keV 
luminosity detection limits of psIO 41,0 and sslO 41 ' 5 ergs s _1 , respectively, at z = 0.7. Such sensitivity 
limits are sufficient to detect moderately-powerful AGNs and normal, X-ray luminous early-type 
galaxies (e.g., NGC 1600; Sivakoff etal. 2004) out to z = 0.7. 

Recently, Mcintosh et al. (2005; hereafter, M05) isolated a sample of 728 optically selected 
early-type galaxies within the E-CDF-S from the GEMS (Galaxy Evolution from Morphologies 
and SEDs; Rix et al. 2004) and COMBO-17 (Classifying Objects by Medium-Band Observations 
in 17 Filters; Wolf etal. 2004) surveys, which overlap with «0.23 deg 2 (»77%) of the E-CDF-S. 
The M05 sample spans a redshift range of z « 0.2-1.0, and due to its depth and relatively large 
comoving volume, it is ideal for investigating the evolution of early-type galaxies over the last half 
of cosmic history. In this paper, we utilize the M05 sample and a supplementary sample of 64 
additional z « 0.1-0.2 early-type galaxies located within the E-CDF-S (selected using the same 
techniques discussed in M05) to characterize the redshift evolution of the X-ray properties of early- 
type galaxies. Our aim is to measure the average X-ray emission, via stacking techniques (applied to 
X-ray detected and X-ray undetected sources), from distant normal early- type galaxy populations. 
We investigate separately the cosmic-time evolution of (1) the hot interstellar gas content through 
optically luminous galaxy samples and (2) the LMXB populations through optically faint galaxy 
samples. 

The Galactic column density is N u ~ 8.8 x 10 19 cm" 2 for the E-CDF-S (Stark etal. 1992). 
All of the X-ray fluxes and luminosities quoted throughout this paper have been corrected for 
Galactic absorption. Unless stated otherwise, we make reference to optical magnitudes using the 
Vega magnitude system. Hq = 70 km s _1 Mpc -1 , Qm = 0.3, and = 0.7 are adopted throughout 
this paper (e.g., Spergel etal. 2003), which imply a look-back time of 6.3 Gyr at z = 0.7. 

2. Data Analysis 

2.1. Sample Selection 

We started with the M05 sample of 728, z = 0.2-1.0 early-type galaxies and a supplemen- 
tary sample of 64 z ~ 0.1-0.2 galaxies within the E-CDF-S field (i.e., 792 total galaxies), which 
were selected following the same methods outlined in M05 using the same data. The total sam- 
ple was selected using (1) rest-frame U — V red-sequence colors derived from COMBO-17 and (2) 
quantitative optical morphologies (via the Sersic indices; n > 2.5) as observed from HST imag- 
ing through GEMS (see § 2.3 of M05 for details). The sample generated using these criteria is 
representative of the early-type galaxy population as a whole and is highly complete down to the 
selection limit of R k, 24. Photometric redshifts for these early-type galaxies were provided by 
COMBO-17 and are used throughout this paper. Using secure redshifts from various spectroscopic 



- 7- 



surveys in the E-CDF-S region (e.g., Le Fevre etal. 2004; Szokoly etal. 2004; Mignoli etal. 2005; 
Vanzella et al. 2005, 2006; Silverman et al. 2006), which were available for «20% of the M05 sample, 
we confirmed that the photometric redshifts for our sample are highly accurate (i.e., ~50% of the 
galaxies have Sz/[1 + z spcc ] < 0.02 and ~75% have Sz/[1 + z spec \ < 0.03). The selection of early- 
type galaxies via these quantitative methods isolates galaxies with visual morphologies ranging 
from E-Sa galaxies, and a large majority of these are E/S0 galaxies (Bell etal. 2004b). We note 
that the X-ray properties of E and SO galaxies are observed to be similar (e.g., Forman etal. 1985; 
Blanton et al. 2001; Xu et al. 2005). 

Our primary goal was to investigate the redshift evolution of the X-ray properties of nor- 
mal early-type galaxies. To achieve this better, we filtered our original sample to include 539 
galaxies within the redshift range z ~ 0.1-0.7. This range was chosen to optimize our use of the 
E-CDF-S X-ray sensitivity to allow for the detection and identification of low-luminosity AGNs 
over a broad range of cosmic look-back times (< 6.3 Gyr). As discussed in § 1.3, the Chandra ob- 
servations of the E-CDF-S are sufficiently sensitive to detect AGNs with 0.5-2.0 keV luminosities 
Lx ~ 10 415 ergs s _1 over this entire redshift range (z ^ 0.7). Furthermore, given that the X-ray 
luminosities of AGN-hosting optically luminous (Lb ~ £_B, C rit) early-type galaxies in the local uni- 
verse span Lx ~ 10 40 - 3-42 - 9 ergs s _1 (mean Lx ~ 10 42 ergs s -1 ; O'Sullivan etal. 2001), this limit 
should adequately prevent powerful AGNs from dominating stacked photon counts (assuming the 
AGN activity is not widespread). 



2.1.1. Removing AGN Contamination 

We utilized multiwavelength observations of the E-CDF-S region to obtain a census of the 
X-ray-detected AGNs and normal galaxies within our sample. We matched galaxies in our sam- 
ple to X-ray sources in the catalogs presented in Alexander etal. (2003) 1 for the ssl Ms CDF-S 
observations and Lehmer et al. (2005b) 2 for the ~250 ks E-CDF-S observations. For a successful 
match, we required that the optical and X-ray centroids be displaced from each other by less than 
1.5 x the radius of the Chandra positional error circles (80-90% confidence), which are given in 
each respective X-ray catalog. We note that investigations of luminous off-nuclear X-ray sources, 
which may lie outside our matching radius, have shown that there is a dearth of such sources within 
early-type galaxies (e.g., Irwin etal. 2004; Lehmer etal. 2006). The Chandra source catalogs were 
generated using wavdetect (Freeman et al. 2002) at false-positive probability thresholds of 1 x 10 -7 
and 1 x 10 -6 for the ~1 Ms CDF-S and the ~250 ks E-CDF-S, respectively. However, as demon- 
strated in § 3.4.2 of Alexander etal. (2003) and § 3.3.2 of Lehmer etal. (2005b), legitimate lower 



1 See http://www.astro.psu.edu/users/niel/hdf/hdf-chandra.html for the relevant source catalogs and data prod- 
ucts for the «1 Ms CDF-S. 

2 See http://www.astro.psu.edu/users/niel/ecdfs/ecdfs-chandra.html for the relevant source catalogs and data 
products for the ^250 ks E-CDF-S. 
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significance X-ray sources, detected by running wavdetect at a false-positive probability threshold 
of 1 x 10~ 5 , can be isolated by matching with relatively bright optical sources. Since the sky surface 
density of z < 0.7, R < 24 early-type galaxies is relatively low (psO.65 galaxies arcmin -2 ), we can 
effectively use this technique to isolate X-ray sources within our sample. We estimate that when 
using wavdetect at a false-positive probability threshold of 1 x 10~ 5 , we expect ss0.5 spurious 
sources. In total, we detected 49 early-type galaxies in at least one of the 0.5-2.0 keV, 2-8 keV, or 
0.5-8.0 keV bandpasses. 

We used the following criteria to identify candidate AGNs within our sample: 

1. Hard X-ray Emission: Our best discriminator of obscured (AT H > 10 22 cm~ 2 ) AGN activity 
is the presence of a hard X-ray spectrum, which is characterized by a relatively shallow 
effective X-ray spectral slope (i.e., T c g 1.5). We classified X-ray sources that were detected 
only in the 2-8 keV bandpass as potential AGN candidates. For sources detected in both 
the 0.5-2.0 keV and 2-8 keV bandpasses, we required that the hardness ratio measured using 
these bandpasses, $2-8 kcv/^0.5-2.0 keV (where $ is the observed count-rate in each respective 
bandpass), be greater than 0.5 (corresponding to an effective photon index of T e fj ^ 1.5); these 
sources would have spectral slopes shallower than those expected for a pure LMXB population 
(e.g., Church k, Balucihska-Church 2001; Irwin etal. 2003). We found that all sources with 
2-8 keV detections were classified as AGN candidates. Furthermore, there were a few sources 
that were detected in only the 0.5-8.0 keV bandpass. We argue that since these sources were 
not detected in the 0.5-2.0 keV bandpass, our most sensitive bandpass, then there must be 
a significant contribution from the 2-8 keV bandpass such that $2-8 kev/^0.5-2.0 keV > 0.5. 
These sources were flagged as AGN candidates and removed from our stacking analyses. Using 
criterion 1, we identified a total of 29 AGN candidates. 

2. X-ray-to-Optical Flux Ratio: Spectral hardness is a good indicator of obscured AGN activity; 
however, powerful unobscured AGNs with steep spectral slopes (r ~ 2) could still be missed 
if this were our only means for identifying candidate AGNs. Another useful discriminator 
of AGN activity, which aids in the identification of luminous unobscured AGN candidates, 
is the 0.5-8.0 keV X-ray-to-optical flux ratio, /0.5-8.0 kev//R ( e -g-> Maccacaro etal. 1988; 
Hornschemeier etal. 2000; Bauer etal. 2004). We use the criterion log(/o.5-8.o kev/fn) > — 1 
(see § 4.1.1 of Bauer etal. 2004 for justification) to identify unobscured AGN candidates in 
our sample; two candidates, which were not identified by criterion 1 above, were found using 
this criterion. We note that 16 out of the 29 AGNs that satisfied criterion 1 (f»55%) also 
satisfied criterion 2. 

3. Radio-to- Optical Flux Ratio: We used radio maps (1.4 GHz; limiting la depth of «14^Jy) 
from the Australia Telescope Compact Array (ATCA; PI: A. Koekemoer; Afonso etal. 2006), 
which cover the entire E-CDF-S region, to identify additional potential AGNs. We matched 
the positions of radio-detected sources to those of our early-type galaxies using a 2'.'0 match- 
ing radius and found 18 matches. Of the 18 radio-detected sources in our sample, we find 
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that three are detected in the X-ray bandpass. Of the three X-ray-detected sources, only 
CXOECDFS J033228.81-274355.6 was not previously identified as an AGN candidate using 
criteria 1 and 2; however, visual inspection of the radio source coincident with CXOECDFS 
J033228. 81-274355. 6 reveals that it has a clear FR II radio morphology. We find that the 
three AGN-like X-ray-detected radio sources have relatively large radio-to-optical flux ratios, 
log(i//„[1.4 GHz] /vf v [5000 A]), as expected for AGNs (e.g., Kinney etal. 1996). Out of the 
remaining 15 radio-detected matches, 13 sources that were not detected in the X-ray band- 
pass were found to have uncharacteristically high radio-to-optical flux ratios for early-type 
galaxies, log(i//„[1.4 GHz] /vf v [5000 A]) > -4. These 13 sources were classified as potential 
AGNs and were excluded from our stacking analyses described below. We note that in ad- 
dition to AGNs, star-forming galaxies (e.g., edge-on dusty galaxies) that were misclassified 
as early- type galaxies may also satisfy this criterion (see, e.g., Fig. 8 of Barger etal. 2006); 
therefore, this crierion also helps to guard against the inclusion of such sources. 

To summarize, we have detected 49 early-type galaxies in the X-ray band and have identified 32 
X-ray-detected AGN candidates: 29 from criterion 1, two additional sources from criterion 2, and 
one additional source from criterion 3. Finally, using criterion 3 we identified an additional 13 AGN 
candidates that were not detected in the X-ray band. Therefore, we identified a total of 45 AGN 
candidates from our sample of 539 early-type galaxies. The remaining 17 X-ray-detected sources 
that we do not classify as AGN candidates are considered to be normal early-type galaxies. These 
normal galaxies are included in our stacking analyses, and as we discuss in § 2.2 below, including 
these galaxies does not significantly affect our results. A more detailed discussion of the X-ray 
properties of our X-ray-detected sources is presented in § 3.1 below. We note that the three criteria 
discussed above may not be completely sufficient to classify all X-ray-detected sources that are truly 
AGNs as AGN candidates (see, e.g., Peterson etal. 2006). Such a misclassification is possible for 
low-luminosity AGNs that are only detected in the more sensitive 0.5-2.0 keV bandpass, which have 
2-8 keV emission too weak for an accurate classification. Furthermore, we also expect there to be 
AGNs that lie below the X-ray detection threshold that are not identified here. However, in § 3.2.2, 
we use the 2-8 keV AGN fraction as a function of X-ray luminosity to argue quantitatively that we 
do not expect misclassified AGNs (detected only in the 0.5-2.0 keV bandpass) and low-luminosity 
AGNs below the X-ray detection limit to have a serious impact on our results. 



2.1.2. Normal Early-Type Galaxy Samples 

A large majority of the normal galaxies in our sample were not detected individually in the deep 
Chandra observations over the E-CDF-S. We therefore implemented stacking analyses (see § 2.2 for 
a description) of galaxy samples, which were constructed to cover two optical luminosity ranges in 
multiple redshift bins. We first divided the original sample of 539 early-type galaxies by luminosity 
into optically luminous {Lb ~ 10 10 ~ n Lb,®) and optically faint {Lb ~ 10 9 ' 3 ~ 10 Lb,®) samples. 
The two luminosity ranges were motivated by the observed physical differences between local early- 
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type galaxies within these two luminosity ranges (see discussion in § 1.1) and optical completeness 
limitations of the M05 sample (i.e., the R < 24 completeness limit). As mentioned in § 1.1 and 
displayed in Figure 1, the optically luminous galaxies have soft X-ray spectra dominated by hot 
interstellar gas, while the optically faint early-type galaxies have X-ray spectra powered mainly by 
LMXBs with a small contribution from hot interstellar gas. Galaxies from both luminosity classes 
have hard X-ray emission dominated by LMXBs. 

We divided the optically luminous sample into four redshift bins of equal comoving volume 
(z « 0.10-0.41, 0.41-0.53, 0.53-0.62, and 0.62-0.70), and due to optical completeness limitations, 
we divided our optically faint sample into two redshift bins (z ps 0.10-0.41 and 0.41-0.53). Al- 
though the E-CDF-S subtends a relatively small solid angle, we note that each redshift bin has a 
radial depth that exceeds 500 Mpc, which is a factor of ^ 5 times larger than the extent of the 
largest structures in the Universe (e.g., Springel etal. 2006). We created histograms of the galaxy 
luminosity-distance distributions on scales of «100 Mpc and did not observe any dominating den- 
sity "spikes." This suggests that our samples are not being significantly affected by cosmic variance. 
Furthermore, provided the physical properties of field early-type galaxies are not strongly affected 
by their environments in a systematic way, we suggest that cosmic variance should not affect our 
results even if present. In Figure 2a, we show rest-frame I?-band luminosity, Lb (in solar units; 
Lb,q = 5.2 x 10 32 ergs s _1 ), versus redshift for the 539 galaxies (dots and crosses) that make 
up our "general sample"; we indicate the divisions of the sample by optical luminosity (shaded 
regions) and redshift (vertical dotted lines). Sources highlighted by large circles (both open and 
filled) are those that are individually detected in X-rays; the open circles indicate X-ray-detected 
AGN candidates, and the filled circles indicate the remaining detected sources, which are classified 
as normal early- type galaxies. 

In addition to the general sample, we have created "faded samples" of early-type galaxies. 
The faded samples were constructed using the general sample discussed above; however, we have 
corrected the rest-frame -B-band luminosities for passive evolution. These corrections take into 
account the fact that the optical power output of early-type galaxies has been fading with cosmic 
time due to stellar evolution (see § 1). For each galaxy, we have computed an evolved, z = 
.B-band luminosity (Lb,o) following the best-fit optical-size-dependent solutions for optical fading 
presented in Tables 1 and 2 of M05. M05 present eight best-fit relations for AMy(z) (i.e., the 
difference between the absolute U-band magnitude at redshifts z and z = 0) versus z, which differ 
in their evolutionary scenarios (based upon formation scenarios and the PEGASE models of Fioc &i 
Rocca-Volmerange 1997) and whether or not the fit is constrained to AMy(z = 0) = 0. Throughout 
this paper we quote results based upon the best-fit solution that assumes AMy(z = 0) = and 
a single-burst evolutionary model with a formation redshift of Zf orm = 3 and a metallicity of 
[Fe/H] = -0.2 (see § 3.2 of M05 for details). We converted the V-band relation to the -B-band 
following AMb ~ 1.1 x AMy. Additional X-ray analyses, using the seven alternative evolutionary 
scenarios presented in M05, were also performed, but no material differences were observed in our 
results. Figure 2b shows Lb,o versus z for our faded sample; symbols, lines, and shaded regions have 
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the same meaning as in Figure 2a. In generating samples for stacking, we made the same sample 
divisions (of both luminosity and redshift) as we did for the general sample, except we used Lb,o 
to discriminate between optically luminous and faint galaxies. This approach allows us to place 
constraints on the evolution of the X-ray activity from distant galaxies selected from the optical 
band as they would appear today. Furthermore, based on local estimates of stellar mass-to-light 
ratios, we can use Lb,o to estimate the stellar masses of these systems. Using a Chabrier initial mass 
function (Chabrier 2003), we estimate that the luminosity ranges Lb,o = 10 9,3-10 and 10 10 ~ n Lb,© 
correspond roughly to stellar-mass ranges of w ^q9.9-io.6 an( ^ ^ jq10.6-ii.6 respectively. In 
the presentation below, we focus our attention on results drawn from the faded samples since these 
are expected to be the most physically meaningful. 



2.2. X-ray Stacking Technique 

In order to address the fact that the X-ray emission from early-type galaxies is dominated 
by different physical processes in different energy bands, we performed stacking analyses in three 
bandpasses: 0.5-1.0 keV (SB1), 3 0.5-2.0 keV (soft band; SB), and 2-8 keV (hard band; HB). 
The shaded bars in the panels of Figure 1 show the spectral coverage of these bandpasses at the 
median redshifts of our optically luminous (^median ~ 0.55) and faint (^median ~ 0.39) samples for 
the X-ray spectral energy distributions (SEDs) of NGC 1600 (top panel) and NGC 4697 (bottom 
panel), respectively. For optically luminous galaxies, the SB1 and SB will effectively sample X-ray 
emission dominated by hot interstellar gas, which produces «80% and «70% of the total emission 
in each respective bandpass. In contrast, 90% of the HB flux originates from LMXB emission. 
For optically faint galaxies, SB and HB will effectively sample the X-ray emission from LMXBs; 
however, the X-ray emission observed in SB1 has roughly equal contributions from hot gas and 
LMXB emission. In our analyses, we used data products from Alexander et al. (2003) for the 
wl Ms CDF-S and Lehmer etal. (2005b) for the «250 ks E-CDF-S (see footnotes 1 and 2). The 
data products are publicly available for all energy bands except for SB1 in the f»250 ks E-CDF-S 
regions, which were generated specifically for the analyses here using the same methods described 
in Lehmer etal. (2005b). The methodology of our stacking procedure was similar to that outlined 
in Lehmer etal. (2005a); for completeness, we describe this procedure below. 

We maximized our stacked signal by optimizing our choices of circular stacking aperture radius 
(from which we extract photon counts for both sources and their backgrounds) and inclusion radius 
(i.e., the maximum off-axis angle within which we include sources for stacking). This process is 
needed because the Chandra point spread function (PSF) increases in size with off-axis angle, 
which degrades the sensitivity for those sources that are far off-axis. For the optimization process, 
we stacked all early-type galaxies in our sample that were not detected individually in the X-ray 
band. In order to obtain a clean X-ray signal, we excluded galaxies that were located ^ 10" from 



3 The 0.5-1.0 keV bandpass was originally denned as SB1 in § 3.1 of Alexander etal. (2003). 
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unrelated sources in the X-ray source catalogs and within the extent of extended X-ray sources, 
which are likely associated with galaxy groups or poor clusters (see § 3.4 of Giacconi et al. 2002 and 
§ 6 of Lehmer etal. 2005b). Sources that lie within both the sal Ms CDF-S and «250 ks E-CDF-S 
were stacked using both observations, as long as the off-axis angle for each observation was within 
our chosen inclusion radius. 

Our optimization procedure was a two-step iterative process that was performed using the 
SB, the bandpass in which our signal was maximized. In the first step of the optimization process 
(step one) we held the inclusion radius (i?i nc i) fixed and stacked sources using a variety of circular 
stacking apertures of constant radii (i.e., we did not vary the aperture radius as a function of 
off-axis angle). We used 15 different circular stacking apertures with radii in the range of / .'5-3'. / 
to obtain a relation for how the signal-to-noise ratio (S/N) varied as a function of aperture size. 
For a given aperture, we stacked the photon counts and effective exposure times from each galaxy 
position and summed them to obtain total source-plus-background counts, S, and a total exposure 
time, T, respectively. We estimated total background counts, B, using Monte Carlo simulations 
and the background maps described in § 4.2 of Alexander etal. (2003) for the «1 Ms CDF-S and 
§ 4 of Lehmer etal. (2005b) for the «250 ks E-CDF-S. We shifted our aperture from each galaxy 
position randomly within a circular region of radiud ~25", extracted background counts from each 
relevant background map, and summed the counts to obtain an estimate of the total background 
counts. This procedure was repeated 10,000 times to obtain an estimate of the local background 
and its dispersion. Our best estimate of the total background counts, B, was approximated by 
taking the mean background calculated from the 10,000 Monte Carlo trials. For each of the 15 
circular apertures, we computed the relevant signal-to-noise ratio (S/N) using the equation, S/N 
= (S — B)/s/B. After performing the stacking procedure for all 15 different apertures, we identified 
the aperture radius that produced the maximal S/N. In the second step of the optimization process 
(step two), we held the optimized aperture determined from step one fixed, but this time we 
stacked sources using 15 different inclusion radii ranging from 1^5—1 lfO to obtain S/N as a function 
of inclusion radius. After performing the stacking procedure for all 15 different inclusion radii, we 
identified the optimal inclusion radius. For a clean stacking signal, we expect that S/N will be 
proportional to the inclusion radius (i.e., S/N oc [Number of sources] 1 / 2 oc J 1 I 2 oc R- nic \), and 

therefore we chose the optimal inclusion radius based on where the S/N-i?i nc i relation appears to 
deviate significantly from linear. This method helps to guard against the inclusion of contaminating 
AGNs in the low-sensitivity regions at large off-axis angles. Using the optimal inclusion radius 
determined with this method, we repeated step one. This two-step process was run iteratively until 
a converged solution was obtained. We found that a stacking aperture of ~1''5 and an inclusion 
radius of ~7'0 produced the optimal signal; we note that this choice of optimized i?i nc i is somewhat 
smaller than that determined in Lehmer et al. (2005a) due to differing X-ray exposures and the 
more conservative approach that we have adopted for choosing the optimal i?i nc i. Hereafter, we 
utilize these values for stacking aperture and inclusion radius in our analyses. In Figure 3, we show 
S/N as a function of inclusion radius (Fig. 3a) and aperture size (Fig. 3b), which were obtained by 
holding the aperture radius constant at l'/5 and the inclusion radius constant at 7f0, respectively. 
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After excluding sources that were (1) classified as AGN candidates (via the three criteria 
outlined in § 2.1.1), (2) located at off-axis angles greater than 7'.0, (3) within regions of extended 
X-ray emission, and (4) within 10" of an unrelated X-ray-detected sources, we were left with 
general and faded samples of 276 and 229 galaxies, respectively; in both samples, 13 galaxies 
were X-ray-detected sources. Using these samples, we stacked the X-ray properties following the 
procedure described above. The sizes of the general and faded samples differ most notably because 
of the fading of high-redshift optically luminous galaxies out of the selected luminosity range. We 
stacked samples both with and without X-ray-detected normal galaxies included, and no material 
differences were observed in our results; therefore, X-ray-detected sources that were classified as 
normal galaxies were included in our stacking analyses. Figure 4 shows the spatial distribution 
of the 229 stacked sources in our faded sample over the E-CDF-S region. Sources stacked using 
the «1 Ms and «250 ks observations are shown as circles and diamonds, respectively; the sources 
indicated with filled circles were stacked using both observations. For each stacked sample, we 
required the S/N be greater than or equal to 3 (i.e., ^ 99.9% confidence) for a detection. For 
stacked samples without significant detections, 3<r upper-limits were placed on the source counts. 

Using net counts (i.e., S — B) from our stacked samples and adopted X-ray SEDs (see below), 
we calculated absorption-corrected fluxes and rest-frame luminosities. Due to the fact that our l'/5 
radius stacking aperture encircles only a fraction of the PSF 4 for sources at relatively large off-axis 
angle, we calculated aperture corrections for each stacked source i. 

Since we are calculating average X-ray counts from the summed emission of many sources of 
differing backgrounds and exposure times, we used a single, representative exposure-time-weighted 
aperture correction, £. This factor, which was determined for each stacked sample, was calculated 
following: 

where Tj is the effective exposure time for each stacked source. The average aperture corrections 
(£) for our samples were «1.5 and wl.8 for the 0.5-2.0 keV and 2-8 keV bands, respectively. We 
note that at the mean redshifts of our samples, the ~1'.'5 radius aperture corresponds to projected 
physical radii in the range of 5.8-10.5 kpc. For early-type galaxies in the local universe, extended 
X-ray emission originating outside of these radii varies considerably among galaxies and generally 
represents only a small fraction (ral-10%) of the total flux (e.g., Fukazawa et al. 2006). Therefore, we 
do not make additional aperture corrections to account for extended X-ray emission. We estimated 
mean observed count rates ($>) using the following equation: 



4 For SB and SB1, the encircled-energy fraction of a l'.'5 radius circular aperture varies from ~100% at off-axis 
angle 9 < 3' to w45% at 9 « 7'. For the HB, this fraction varies from ^80% at 6 < 3' to f«40% at 9 w 7'. 
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* = £{S-B)/T, (2) 

where S, B, and T are defined above. To convert count rate to flux, we used the Galactic column 
density given in § 1.3 and the best-fit X-ray spectral energy distribution (SED) for NGC 1600 for 
the optically luminous galaxies (see Fig. la) and a power-law SED with T = 1.8 for the optically 
faint galaxies (see § 1.1 for justification). We note that for a range of reasonable SED choices, we 
find systematic fractional uncertainties in the count-rate to flux conversion of ^ 10% for SB and 
HB and «50% for SB1. Mean rest-frame X-ray luminosities Lx,R were calculated following: 

<Lx,R)~47r(4)</x,o>#, (3) 

where di is the luminosity distance, (/x,o) is the mean observed-frame X-ray flux, and K is a 
unitless conversion factor, which relates the observed-frame X-ray flux to the rest-frame luminosity 
using our adopted SEDs. The fractional errors on (Lx,r) due to uncertainties in {d 2 L ) range from 
«10% to wl% for our optically luminous z « 0.10-0.41 and z 0.62-0.70 samples, respectively. 
Due to the relatively broad energy ranges our bandpasses encompass and the relatively small 
spectral shifts observed over the redshifts of sources in our samples, we used observed-frame fluxes 
to compute rest-frame luminosities. For our adopted SEDs, we found that values of K varied 
between 0.8 and 1.1, depending on the bandpass and redshift of the source. For all bandpasses, 
SED choice contributes small systematic uncertainties in K, which range from « 10-20% depending 
on redshift. We note that the mean X-ray luminosity is expected to be closely representative of 
a typical galaxy within the confined optical luminosity ranges used here. Moreover, using the 
O'Sullivan etal. (2001) sample of local early-type galaxies, we find that L£ can /L x lcdian ~ 1.7 an d 
1.2 for optically luminous and faint samples respectively. Since we are only able to calculate mean 
quantities via the X-ray stacking used here, we utilize mean quantities throughout. Results from 
our stacking analyses are presented below in § 3.2. 



3. Results 

3.1. Individually Detected X-ray Sources 

In Table 1, we present the properties of the X-ray-detected early- type galaxies. X-ray source 
properties were determined following the methods outlined in Alexander et al. (2003) and Lehmer 
et al. (2005b) for the CDF-S and E-CDF-S, respectively. Using the matching criterion discussed 
in § 2.1.1, we matched 39 of our early- type galaxies to sources that were included in the main 
Chandra catalogs of either Alexander etal. (2003) or Lehmer etal. (2005b), and an additional ten 
matches were identified using X-ray sources detected using wavdetect at a false-positive probability 
threshold of 1 x 10~ 5 (i.e., 49 total detected sources). Figure 5 shows the i?-band magnitude (from 
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COMBO-17) versus 0.5-8.0 keV flux for all 49 X-ray detected early-type galaxies in our sample; nor- 
mal galaxies and AGN candidates are indicated as filled and open circles, respectively. As discussed 
in § 2.1.1, we classified 32 («65%) X-ray-detected sources as AGN candidates. We note that the 
majority of the AGN candidates have log /o. 5-8.0 kev//R > — 1> however, a few of these candidates 
have log ib.5-8.0 kcv//R < -2, including the FR II source CXOECDFS J033228.81-274355.6. AGNs 
with log /o. 5_8.o kcv//R < —2 are either significantly obscured or relatively X-ray- weak AGNs. 

As discussed in § 1.2, transient AGN feedback may play an important role in heating the hot gas 
in normal early- type galaxies. In order to understand and constrain the influence of transient AGN 
activity within early-type galaxies, we computed the X-ray-luminosity-dependent cumulative AGN 
fraction, fc (i.e., the fraction of early- type galaxies harboring an AGN with a 2-8 keV luminosity 
of L2-8 keV or greater). We made use of the 2-8 keV bandpass because of its ability to probe 
relatively unattenuated X-ray emission in a regime of the X-ray spectrum where we expect there 
to be minimal emission from normal galaxies (see also § 2.1.1 for more details). Figure 2 illustrates 
that the majority of the AGNs in our samples (open circles) originate within optically luminous 
(Lb ~ 10 10 Lb,q) early-type galaxies, and therefore when computing fc, we used only optically 
luminous galaxies; we note that the number of AGNs within our optically faint samples is too low 
to obtain a respectable constraint on fc- In order to quantify the redshift evolution of fc, we 
split our optically luminous samples into two redshift intervals of roughly equal size (z ~ 0.10-0.55 
[^median = 0.42] and z Tti 0.55-0.70 [^median = 0.65]). Only two redshift intervals were chosen due to 
statistical limitations on the number of detected AGNs. We computed fc by taking the number of 
candidate AGNs with a 2-8 keV luminosity of L2-8 keV or greater and dividing it by the number of 
early-type galaxies in which we could have detected an AGN with luminosity L2-8 keV- The latter 
number was computed by considering the redshift of each galaxy and its corresponding sensitivity 
limit, as obtained from spatially dependent sensitivity maps (see § 4.2 of Alexander et al. 2003 
and § 4 of Lehmer etal. 2005b); these sensitivity maps were calibrated empirically using sources 
detected by wavdetect at a false-positive probability threshold of 1 x 10~ 5 . Figure 6 shows fc as 
a function of L2-8 keV f° r the two redshift bins considered here. We find suggestive evidence for 
evolution in fc between z ~ 0.42 and z ~ 0.65, which is consistent with the global trend observed 
for luminous AGNs in general (e.g., Brandt & Hasinger 2005). At L2-8 keV ~ 10 42 ' 2 ergs s _1 , where 
fc is most tightly constrained, we find that fc(z = 0.65) = p.l^Vl x fc( z = 0.42). We note 
that although this value is poorly constrained, it is consistent with the (1 + z) 3 evolution observed 
for X-ray luminosity functions of X-ray-selected AGNs (e.g., Ueda et al. 2003; Barger etal. 2005; 
Hasinger etal. 2005), and is in agreement with the stacked constraints on optically-selected early- 
type galaxies set by Brand etal. (2005). We return to the discussion of X-ray-detected AGNs in 
§ 3.2.2 when discussing the undetected AGN contribution to our stacked samples and in § 4.1 when 
discussing the transient AGN contribution to heating the hot interstellar gas in optically luminous 
early-type galaxies. 
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3.2. X-ray Stacking Results 

3.2.1. Stacked Properties 

In Tables 2 and 3, we summarize the average properties of our stacked samples of normal 
early-type galaxies. For illustrative purposes, we created Figure 7, which shows 0.5-2.0 keV (SB) 
adaptively-smoothed stacked images of our faded samples. We detect the average X-ray emission 
from all of our samples in SB, several in SB1, and only two in HB. The two samples from which we 
detect HB emission (the z as 0.65 optically luminous general sample and the z ~ 0.58 optically lu- 
minous faded sample) have HB-to-SBl and HB-to-SB count-rate ratios (i.e., <3?2-8 kev/^0.5-1.0 kcV 
and $2-8 kcv/^0.5-2.0 kev) that are broadly consistent with our adopted X-ray SED. Furthermore, 
all of our stacked samples have mean X-ray-to-optical flux ratios (/o. 5-8.0 kev/f~R.) that are con- 
sistent with those expected for normal galaxies (see col. [5] of Table 3 and Fig. 5). In Figure 5, we 
have plotted mean quantities from our samples as large filled squares and triangles, which represent 
our optically luminous and faint faded samples, respectively; these symbols have been shaded with 
varying grayscale levels to indicate redshift, such that darker shading represents higher redshift 
samples. 

For the purpose of comparing our results with local early-type galaxies, we utilize the O'Sullivan 
etal. (2001; hereafter OS01) sample. The OS01 sample was selected from the Lyon-Meudon Extra- 
galactic Data Archive (LEDA) using morphological type (T < —1.5; E-S0 Hubble types), distance 
(V < 9000 km s _1 ), and apparent magnitude (Bt < 13.5). The LEDA catalog is known to be ~90% 
complete down to Bt = 14. X-ray observations of these galaxies were available mainly from the 
ROSAT PSPC with a significant minority of the data originating from the Einstein IPC. We also 
utilized Chandra data from David et al. (2006) for six of the OS01 galaxies having only X-ray upper 
limits. Figure 8 shows the 0.5-2.0 keV luminosity versus i?-band luminosity for galaxies included 
in the OS01 sample with D < 70 kpc; X-ray luminosities have been normalized to the 0.5-2.0 keV 
bandpass using the X-ray SED adopted in OS01 (i.e., a MEKAL plasma SED with solar metallicity 
and a plasma temperature of 1 keV). We have denoted field, central-cluster, brightest-group, and 
AGN-hosting early-type galaxies as identified from the OS01 sample; several well-studied examples 
have been highlighted (M32, M87, NGC 1399, NGC 1600, NGC 4697, and NGC 5102) for reference. 

As noted in § 1.1, the X-ray and B-band luminosities of local early- type galaxies are observed 
to be correlated, and these correlations follow a power law: 

log L x = a log L B + j3, (4) 

where a and (5 are fitting constants. Using the OS01 sample, we performed linear-regression analyses 
for galaxies with Lb ~ 10 10 Lb,q and Lb ~ 10 10 Lb,q separately to determine a and {3 for each 
luminosity regime. When doing these calculations, we excluded (1) sources at D > 70 kpc, (2) 
sources with X-ray emission that may be significantly influenced by X-ray-emitting gas associated 
with galaxy clusters or groups such as central-cluster and brightest-group galaxies, (3) AGN-hosting 
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galaxies, (4) NGC 5102, due to its anomalously low X-ray luminosity and evidently recent star- 
formation activity (e.g., OS01; Kraft etal. 2005; see Fig. 8), and (5) NGC 4782, which has an 
anomalously large Lb that drives the correlation (i.e., Lb = 10 114 L_b,©; see, e.g., OS01). We 
utilized Kendall's tau (Kendall 1938) to measure the correlation strengths and Buckley-James 
regression (Buckley & James 1979; Isobe et al. 1986) to calculate the best-fit correlation parameters 
(i.e., a and (3) for each luminosity regime. These tools were available through the Astronomy 
SURVival Analysis software package (ASURV Rev. 1.2; Isobe & Feigelson 1990; LaValley etal. 1992). 
We found correlation significances of «6.8<7 and «2.8<t for the Lb ~ 10 10 Lb,q and Lb ~ 10 10 Lb,q 
galaxy samples, respectively. We found that for the optically luminous sample a = 2.61 ± 0.61 and 
(3 = 12.77 (solid line in Fig. 8), and for the optically faint sample a = 1.05 ± 0.27 and (3 = 29.36 
(dashed line in Fig. 8). We note that the Lx~Lb relation for optically faint galaxies is poorly 
constrained with Ri75% of the galaxies having only X-ray upper limits; however, since these galaxies 
are thought to be dominated by LMXBs, investigations of the LMXB luminosity per unit i?-band 
luminosity can offer a consistency check for this relation. Using Chandra observations of 14 E/S0 
galaxies, Kim & Fabbiano (2004) found that for LMXBs with L x > 10 37 ergs s" 1 , (3 = 29.5 ± 0.25 
using a fixed slope of a = 1. This value is the most tightly constrained Lx(LMXB)-L# relation 
available at present and is consistent with other investigations of the discrete-source contribution 
(e.g., Sarazin 1997; OS01 and references therein; Gilfanov 2004) and our calculated Lx~Lb relation 
for local optically faint early-type galaxies. For reference, the expected discrete-source contribution 
from LMXBs is presented in Figure 8 as a shaded region, which represents the dispersion of the 
relation. We note that the Lx~Lb relation for local optically faint early-type galaxies (dashed line) 
is consistent with that expected for X-ray emission originating strictly from LMXBs. 

In Figure 8, we have plotted mean quantities from our samples using the same symbols and 
symbol-shading schemes that were adopted in Figure 5. The plotted error bars for our stacked 
samples were computed by propagating (1) Poisson errors on the source counts (Gehrels 1986), 
(2) 1 a errors on (d 2 L ), and (3) systematic errors on the SED-dependent parameters (i.e., count- 
rate to flux conversion and K); these errors were propagated following the "numerical method" 
described in §1.7.3 of Lyons (1991). In an initial evaluation of these results, we find that the 
average X-ray properties of our optically luminous samples appear to follow the local relation at 
all redshifts. In contrast, our optically faint samples deviate from the local relation significantly 
(1.5cr at z « 0.46) over the redshift range z ~ 0.1-0.5, suggesting there may be some evolution in 
the LMXB populations within these galaxies. 

Figures 9a and 9b show the mean 0.5-2.0 keV luminosities (Lx) of our stacked general and 
faded samples, respectively; optically luminous and faint samples are indicated as dark filled squares 
and triangles, respectively. X-ray-detected sources are shown as circles, which represent both nor- 
mal galaxies (filled circles) and AGN candidates (open circles). We have included the correspond- 
ing mean X-ray luminosities for local early-type galaxies from the OS01 optically luminous and 
faint samples in matched optical-luminosity ranges; as before, these samples were filtered to ex- 
clude galaxies classified as AGNs, central-cluster galaxies, brightest-group galaxies, NGC 5102, and 
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NGC 4782. Mean optical and X-ray luminosities for the OS01 samples were computed using the 
Kaplan-Meier estimator (e.g., Feigelson & Nelson 1986) within ASURV, which appropriately handles 
censored data. When calculating these mean luminosities, we filtered the OS01 sample to include 
only sources with distances < 40 Mpc and < 20 Mpc for the optically luminous and faint samples, 
respectively; these distances represent approximate completeness limits for the optical-luminosity 
ranges used here. We calculated mean X-ray luminosities of log(Lx) ~ 40.2 ± 0.2 and «39.2 ± 0.2 
for the optically luminous and faint samples, respectively; error bars here represent la errors on the 
means (computed with ASURV). These calculations were made using 102 optically luminous galax- 
ies (60 upper limits) and 48 optically faint galaxies (36 upper limits). We attempted to improve 
the X-ray-detection fractions of these samples by using subsamples of galaxies that were created 
with distance limits smaller than the completeness limits quoted above; however, we found no 
improvement in these fractions. 

Figures 10a and 10b show the X-ray-to-S-band mean luminosity ratio (Lx/Lb) for our general 
and faded samples, respectively (symbols have the same meaning as they did in Fig. 9); the expected 
local discrete-source contribution and its dispersion are shown as a horizontal dotted line and the 
surrounding shaded region, respectively. As observed in Figure 8, we find little evidence for evolu- 
tion in our optically luminous samples (both general and faded). Using the available data, we find 
that (Lx/Lb) z =o.7 = [1.0±0.5] and [1.7±0.8] x (Lx/Lb) z =o for our optically luminous general and 
faded samples, respectively. In order to constrain further the allowed redshift evolution of Lx/Lb, 
we utilized the x 2 statistic and a simple evolutionary model, (Lx/Lb) z = (1 + z) v (Lxl 'Lb) z =o- F° r 
this single-parameter model, we constrained rj using 90% confidence errors (Ax 2 = 2.7). We found 
best-fit parameters of ij = — OA^q® (x 2 = 3.4 for 4 degrees of freedom) and r/ = 0.4^q^ (x 2 = 11-3 
for 4 degrees of freedom) for our optically luminous general and faded samples, respectively. For 
the optically faint early-type samples, we observe suggestive redshift evolution in Lx/Lb, and by 
z ~ 0.5 it has increased above the local relation by a factor of ~5.3±4.1 and ~5.6±4.1 for our gen- 
eral and faded samples, respectively. Using the x 2 statistic and the same simple model for redshift 
evolution presented above, we find best-fit values of rj = 4.4 ± 1.7 (x 2 = 1.1 x 10~ 3 for 2 degrees 
of freedom) and r] = 5.4 ± 1.5 (x 2 = 1-0 for 2 degrees of freedom) for our optically faint general 
and faded samples, respectively. We note that the evolution observed for optically faint early-type 
galaxies is largely driven by the value of Lx/Lb at z = 0. Due to the fact that the z = value for 
Lx/Lb is based on 48 sources, with 36 (75%) having only X-ray upper limits, we cannot rule out 
the possibility that Lx/Lb at z = is significantly affected by systematic errors. Moreover, the 
total X-ray emission from optically faint galaxies is expected to vary significantly between galaxies 
due to low numbers of LMXBs and variable amounts of hot gas, and therefore large fractional errors 
are expected for Lx/Lb- We therefore consider this result to be only marginal at present. 

As noted above, the local relation between Lx and Lb is nonlinear for optically luminous 
early- type galaxies. In order to investigate whether such nonlinearities have an effect on our overall 
results, we created Figure 11, which illustrates residuals to the local best-fit relations, logLx = 
a log Lb + 0, for optically luminous (a = 2.61; (3 = 12.77) and faint (a = 1.05; (3 = 29.36) 
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samples. Figure 11 shows that the nonlinearities observed in the local Lx~Lb relation do not 
affect our conclusions drawn from using Lx/Lb as a proxy for evolution. Furthermore, due to 
the approximate equality of the mean values of Lb for all samples of a given luminosity class, 
quantitative analyses that account for nonlinearities in the Lx~L b relations yield roughly identical 
results to those quoted using simply Lx/Lb- 



In this section we assess whether the stacked properties presented above suffer from contamina- 
tion by AGNs with 2-8 keV luminosities below our detection threshold. In § 3.1 (see also Fig. 6), we 
presented the cumulative AGN fraction, fc, the fraction of early-type galaxies harboring an AGN 
with a 2-8 keV luminosity of L 2 -$ kcV or greater. To first order, we can use the functional form of 
fc to generate a census of the AGN population that we expect to be missing due to sensitivity limi- 
tations. As noted in § 3.1, there is evidence that fc evolves with redshift. We modelled this redshift 
evolution of fc using the functional form, fc{z) = (1 + z) 3 fc{ z = 0) (see § 3.1 for justification). We 
also assumed that the dependence of fc upon L2-8 keV was similar at all redshifts over the range 
z ~ 0.1-0.7. Using our faded sample of optically luminous early-type galaxies we computed (fc)z, 
the average cumulative AGN fraction over the redshift range z « 0.1-0.7 (z mc dan = 0.55), following 
the procedure outlined in § 3.1; this was done to obtain a better understanding of the overall shape 
of the fc(L2—8 kcv) curve. Figure 12a (filled circles with error bars) shows our computed values of 
(fc)z as a function of logL 2 _s keV- We fit the (fc)z data points using a quadratic relation (i.e., 
log(/c) 2 = o + ailogL 2 _ 8 kcv + a2(logL 2 _ 8 kcv) 2 , where a = -136, ai = 7.1, and a 2 = -0.1; 
thick solid curve in Fig. 12a) over the luminosity range logi^-s keV = 40-44; this covers the same 
luminosity range for AGNs as the OS01 local sample. Using our best-fit relation for {fc)z an d our 
model for the redshift evolution of fc, we calculated fc(z, log L2-8 kcv) for each of our optically 
luminous faded samples following: 



where z is the median redshift of each sample, z me dian = 0.55 is the median redshift of our best- 
fit redshift-averaged relation, (fc)z- In Figure 12a, we show our estimates of fc for each of our 
optically luminous faded samples; these curves are annotated on the figure. For comparison, we 
have shown the AGN fraction measured for DRGs (see § 1) at z ~ 2.5 by Rubin et al. (2004; open 
diamond) and have extrapolated our model for fc out to z ~ 2.5 (dot-dashed curve). We note that 
even at z « 2.5 our model agrees reasonably well with observations. 

In order to estimate the amount of AGN contamination that may be contributing to our 
stacked signals, it is desirable to convert the cumulative AGN fractions to differential forms (i.e., 
the fraction of galaxies harboring AGNs within discrete X-ray luminosity bins). Using the model 



3.2.2. Assessing Remaining AGN Contamination 
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for fc presented above, we estimated the fractions of galaxies harboring an AGN within luminosity 
bins of width AlogL 2 „ 8 keV = 0-5 for each of our optically luminous faded samples; we refer 
to these as differential AGN fractions, fr>. The histograms in Figure 12b show our estimates of 
fo- Due to the deep Chandra coverage in the E-CDF-S region, a large fraction of the luminous 
AGNs (L 2 -8 keV ~ 10 41,5 ergs s" 1 ) would have been removed from our samples before stacking (see 
§ 2.1.1). We are therefore only interested in the fraction of galaxies falling below our sensitivity 
limit. Using the sensitivity maps described in § 3.1, we determined the fraction of optically luminous 
galaxies within each of our stacked faded samples for which we could not have detected an AGN 
of luminosity keV if present; we refer to these fractions as fg (i.e., the fractions of galaxies 

below our sensitivity limit) and show them in Figure 12c. For each of our optically luminous 
faded samples we calculated the fraction of sources that harbor an undetected AGN with 2-8 keV 
luminosity L 2 _ 8 keV ( in bins of width AlogL 2 -8 keV = 0.5), fa, by multiplying f D by f B . In 
Figure 12d, we show fjj as a function of logL 2 _8 keV- Using fu, we computed the approximate 
AGN contamination with the following summation: 



where the summation is over all bins of AlogL 2 _s keV = 0.5. We find L 2 _§ k c v (contain) rj 
-Lq39. 5-40.3 er g g g -i £ or our sanr ples. F or our samples with z £ 0.45, AGNs with L 2 _s keV ~ 10 41 ergs s 
contribute ^ 70% of the total L 2 _s keV (contam) estimate. Furthermore, when extrapolating our 
model for fc down to much lower values of L 2 _ 8 k e v> we find no significant difference in our esti- 
mates of L 2 s keV (contam). We note that this result is mildly dependent on our extrapolation of fc 
to values of L 2 _ 8 keV ~ 10 39 ~ 40 - 5 ergs s~ 4 . If there exists a large population of radiatively-inefficient 
low-luminosity AGNs in early- type galaxies that radiate within this X-ray luminosity range (e.g., 
advection-dominated accretion flows [ADAFs]), then we may be underestimating L 2 _g keV (contam). 
However, Chandra observations of ~50 early-type galaxies in the local universe have revealed that 
the majority of the central supermassive black holes in these galaxies are radiating at extremely 
low efficiencies and typically have observed Lx ~ 10 39 ergs s" 1 (e.g., Loewenstein et al. 2001; David 
etal. 2005; Pellegrini 2005). Therefore a change in the shape of fc at L 2 _g keV ~ 10 39-40 ' 5 ergs s~ 4 
is not expected. 

We converted L 2 s kev( con tam) to a 0.5-2.0 keV luminosity assuming a power-law model 
with an effective photon index T e g, which was determined by stacking all 24 AGN candidates 
in our optically luminous faded sample; in this calculation we purposely made no attempt to 
correct for intrinsic absorption. For these 24 AGN candidates, we find that the mean band ratio 
^2-8 kcv/^o. 5-2.0 keV = 1-12 ± 0.09, which corresponds to an effective photon index of r e ff = 
0.77 ± 0.08; this effective photon index is in agreement with the effective photon indices measured 
for sources in the Chandra Deep Field-North (CDF-N; see Fig. 14 of Alexander et al. 2003) with 
similar fluxes. If we assume F e g = 0.8 is a characteristic effective photon index for the AGNs we 
expect to be missing, then we find that AGN contamination can account for «4% and pss11% of the 
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observed 0.5-2.0 keV emission from our optically luminous faded samples at z ~ 0.25 and z pa 0.66, 
respectively; this amount does not significantly affect the results presented in § 3.2. f above. We 
note that the X-ray SED used in this calculation has an important effect on the overall estimate 
of the AGN contamination. Since our estimate for contamination decreases as r e g decreases, the 
amount of contamination in our samples may be affected if our choice of r e g- is too flat. However, 
if we choose a steeper effective photon index such as r e fr = 1.4 (the observed spectrum of the X-ray 
background), we still find that the estimated AGN contribution to our 0.5-2.0 keV signal is too low 
( 25% at z sa 0.66) to make a substantive difference to our results. Furthermore, as discussed in 
§ 2.1.1, we have taken additional precautions to eliminate several AGN candidates that were not 
detected in the 2-8 keV bandpass; these sources would not be taken into account in this estimate 
for AGN contamination. We also note that similar results are found when performing the above 
analyses using our optically luminous general samples. 

Due to poor statistical constraints on AGN activity in optically faint early-type galaxies, an 
estimate of fc could not be determined reliably using the present data. Studies of the AGN fraction 
as a function of galactic stellar mass have found that AGNs are much more commonly observed in 
massive galaxies than lower-mass galaxies (e.g., Kauffmann et al. 2003). We therefore approximated 
a strict upper limit to the AGN contribution to the X-ray emission from our optically faint early- 
type samples by using the same model for fc presented above for the more massive optically 
luminous galaxies. Using this model, we estimate that AGNs contribute < 20% and < 40% of 
the X-ray emission from our optically faint samples (both general and faded) assuming T c g = 0.8 
and r e ff = 1.4, respectively. These limits slightly reduce the significance of the quoted evolution 
for our optically faint faded samples such that (Lx/Lb) z =o.5 is estimated to be [4.5 ± 3.3] and 
[3.4 ± 2.5] x (Lx/Lb) z =o for r e g- = 0.8 and F e s = 1.4, respectively; however, we note that these 
limits should be regarded as very conservative. 

4. Discussion 

The above results suggest differing evolutionary histories for optically luminous and faint early- 
type galaxies. As shown in § 3.2.2, our results are not expected to be significantly affected by an 
undetected population of AGNs. Therefore, changes in the X-ray emission with redshift are likely 
the result of global changes in the emission from hot interstellar gas and/or LMXBs. In the sections 
below, we discuss possible interpretations of the results for our optically luminous and faint early- 
type galaxies in turn. 

4.1. Optically Luminous Early- Type Galaxies 

As discussed in § 2.2, we chose the SB for our stacking analyses to sample directly the X-ray 
emission from hot interstellar gas. Therefore, the near constancy of Lx/ Lb with redshift can be 
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largely explained as a general balance between the energy losses from the hot gas (A22 gas ) and 
the energy deposition from heating mechanisms (AE^ating) over each cooling time, t coo \. Here we 
investigate the relative contributions from feedback mechanisms to constrain physical models of 
the heating of the hot gas in early- type galaxies. As noted in § 1.2, the typical inferred radiative 
cooling time for the central regions of an optically luminous early-type galaxy is i coo i « 10 s yr. 
Since each of our optically luminous early-type galaxy redshift bins are larger than the cooling 
timescale (i.e., our redshift bins have temporal widths in the range of ps0.5-4.4x10 9 yr), we can 
estimate the redshift-dependent energy components following: 



Here, L gas is the redshift-dependent globally-averaged power output from the gaseous component of 
our optically luminous early-type galaxies. The first two terms of equation 7b represent AGN heat- 
ing from both radiative and mechanical feedback power, respectively. The radiative feedback power 
is represented as the product of the average 2-8 keV AGN luminosity per galaxy (L 2 _ 8 keV,AGN)> 
its bolometric correction factor (7bc ~ 30; e.g., Marconi et al. 2004; Barger et al. 2005), and the effi- 
ciency factor describing the coupling between radiation and the hot interstellar gas (e ra( j). We note 
that in cases where the Compton temperature of the AGN SED falls below the temperature of the 
hot interstellar gas, radiation from the central AGN may effectively cool the gas and thereby drive 
e ra( j to negative values (see, e.g., § 6 of Nulsen & Fabian 2000; Ciotti & Ostriker 2001). Mechanical 
power (e.g., through AGN jets) is likely to be a very important feedback mechanism, and we have 
indicated its contribution as L m cch,AGN- Finally, Aether represents additional energy input from 
alternative forms of heating over each cooling time (see below) . 

Using the hot-gas component of the X-ray SED for NGC 1600 (see dotted curve in Fig. la), we 
estimate L gas ~3x Lo. 5-2.0 kcVi f° r ° ur optically luminous samples, this amounts to a mean value 
of (L gas ) ~ 8 x 10 40 ergs s _1 . We measured L2-8 kcV,AGN using our redshift-dependent model for 
the differential fraction, fjj, which was presented in § 3.2.2 (see also Fig. 12b). Duty cycles for AGN 
activity in any given galaxy are expected to be shorter than the timescales represented by each of our 
redshift bins; however, since we are considering large populations of early-type galaxies, we do not 
expect significant variations in the AGN fraction (as measured from the E-CDF-S "snapshot" ) at 
any given time within each redshift bin. Therefore using this model, we can calculate L 2 -8 kcV,AGN 
following: 
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e r ad7BC-^2-8 keV,AGN*cool + -^mech,AGNi C ool + AE ot her- 



(7a) 
(7b) 




(8) 



1 



where the summation is over bins of A log £2-8 keV = 0.5 and covers the luminosity range L2-8 keV = 
10 40-44 ergs s _1 . We found L2-8 kcV.AGN ~ 9.1 and «21 xlO 40 ergs s" 1 per galaxy at z ~ 0.25 
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and z 7n 0.66, respectively. In Figure 13, we show L 2 _8 keV,AGN as a function of redshift for 
our optically luminous faded samples (filled squares). For comparison, we also show stacking 
results from the Brand etal. (2005) samples of z ~ 0.4-0.9 early-type galaxies (open diamonds), 
which have mean i?-band absolute magnitudes that are well matched to those of our optically 
luminous faded samples; these mean luminosities are dominated by X-ray emission from AGNs 
and therefore provide a good estimate of L 2 _8 kcV,AGN (~80-90% of the emission is from AGNs). 
Using these data and our model for the evolution of the X-ray emission from transient AGNs, 
L2-8 keV.AGN = L 2 - 8 keV,AGN, z =o (1 + zf , we found L 2 - 8 kcV,AGN,z=o ~5x 10 40 ergs s" 1 (see the 
dotted curve in Fig. 13). 

Using the above relations and the assumption that A_E gas = A-E^eating) we arrive at the fol- 
lowing relation: 

igas = erac«7BC^2-8 keV,AGN,z=0 (1 + z) 3 + £ m ech,AGN + A.E other /i cool (9) 

As discussed in § 1.2, transient AGN feedback is expected to play a significant role in the heating 
of the hot interstellar gas. If we assume that AGN feedback is largely responsible for keeping 
the gas hot, then we can neglect the last term of equation 9 (i.e., AE^thcrAcooi <C L gas ). With 
this assumption and the observed constancy of L gas with redshift, we infer that the strongly- 
evolving radiative power must be poorly coupled to the hot interstellar gas, such that e ra( j <C 0.05. 
This suggests that mechanical AGN power (i.e., -L mec h) dominates the feedback over the redshift 
0.0-0.7 and does not evolve in the same way as the radiative power. We note that 
it is also possible that A.E ot h er may play some non-negligible role in the heating of the hot gas. 
Additional heating sources may include inward thermal conduction from the large reservoirs of hot 
gas found in the outer regions of early-type galaxies (e.g., Narayan & Medvedev 2001; Brighenti & 
Mathews 2003), Type la supernovae and stellar winds (e.g., Loewenstein & Mathews 1987, 1991), 
and infalling circumgalactic gas (e.g., Brighenti & Mathews 1998); however, the influence of these 
heating sources is presently not well constrained. 



4.2. Optically Faint Early- Type Galaxies 

For our optically faint early-type galaxy samples, we found suggestive evidence for redshift 
evolution in Lx/Lb over the redshift range z ~ 0.0-0.5 (see § 3.2.1 and Figs. 10 and 11). Although 
some of the observed emission may be due to AGN activity (§ 3.2.2), there remains suggestive 
evidence that normal X-ray activity is evolving with redshift, and we discuss possible scenarios 
explaining this evolution below. 

One possible driver of Lx/Lb evolution may come from global changes in the LMXB popu- 
lations within optically faint early-type galaxies. As discussed in § 1.2, LMXBs from primordial 
binaries within the galactic field are expected to dominate the overall LMXB emission from optically 
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faint early-type galaxies; this differs from the LMXB emission from optically luminous early-type 
galaxies, which originates primarily from globular clusters. LMXBs from primordial binaries emerge 
in the wake of star-formation epochs ~1-10 Gyr following a major star-formation event. Therefore, 
changes in the mean early-type galaxy stellar age with redshift should result in observed changes 
in the mean X-ray emission from these systems. Furthermore, these changes are expected to be 
most conspicuous «1 Gyr after major star-formation events (e.g., White & Ghosh 1998; Ghosh & 
White 2001). Galaxy formation scenarios that favor a more recent emergence of the optically faint 
early-type galaxy population onto the red-sequence such as downsizing or mass-dependent merging 
histories (see discussion and references in § 1) would predict significant evolution of the LMXB 
emission from these systems. 

A second source of evolution of the X-ray emission from optically faint early-type galaxies 
could in principle come from cooling of the hot X-ray-emitting gas. However, studies of optically 
faint early-type galaxies in the local universe have shown that the hot gas emission generally makes 
up a minority fraction (typically «40%) of the total X-ray emission (e.g., David etal. 2006) and is 
therefore less likely to be completely responsible for the observed evolution than LMXBs. 

4.3. Future Work 

Understanding the evolution of the X-ray properties of early-type galaxies could be greatly 
improved by (1) constraining better the X-ray properties of local optically faint early-type galaxies 
and making a census of their AGN populations, (2) performing additional investigations using 
other deep Chandra fields that have complementary HST coverage, and (3) conducting deeper 
observations using Chandra or future X-ray missions. These possibilities are discussed in more 
detail below. 

In Figure 8, we showed the values of Lx and Lb for local early-type galaxies from the OS01 
sample, which is the largest, uniformly-selected sample available for studying the Lx^Lb correla- 
tions of local early-type galaxies. The majority («75%) of the isolated optically faint early-type 
galaxies have only X-ray upper limits, which has restricted our interpretation of the redshift evo- 
lution of these galaxies. New Chandra observations of the galaxies having only X-ray upper limits 
could not only improve the characterization of the Lx~Lb correlation at lower Lb, but could also 
provide useful insight into the role of hot interstellar gas and LMXBs within the galactic field, 
which are presently not well constrained. Furthermore, stacking analyses of well-chosen samples 
of these galaxies could provide useful statistical insight into the mean X-ray properties of these 
galaxies and mitigate the effects of poor source statistics for individual galaxies. 

Chandra stacking analyses using additional samples of early-type galaxies could effectively 
reduce the sizes of the errors on mean quantities and/or allow for the analyses of samples in 
more finely partitioned bins of redshift and optical luminosity. An important requirement for such 
studies is to obtain an adequate census of the underlying AGN population, which may significantly 
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influence the stacking results. In order to remove effectively AGNs from stacked signals at z ^ 0.7, 
relatively deep Chandra observations are necessary. Using the luminosity-dependent AGN fractions 
determined in § 3.2.2, we suggest that AGNs will provide significant contamination for Chandra 
exposures of ^ 100 ks. Therefore, studies of distant early-type galaxies using multiwavelength data 
(including HST coverage) from already existing Chandra fields such as the ~2 Ms CDF-N or the 
«200 ks All-wavelength Extended Groth Strip International Survey (AEGIS; Nandra etal. 2005; 
Davis et al. 2006) would improve the present situation. 

Finally, deeper Chandra observations of already existing fields (most notably the CDF-N) 
would provide an improved census of the low-to-moderate luminosity AGN population at higher 
redshifts and connect the X-ray properties of these relatively passive early-type populations with 
those of their higher redshift progenitors (e.g., DRGs, EROs, and distant submillimeter-emitting 
galaxies; e.g., Alexander etal. 2005). In addition to the improvement that additional Chandra 
observations could provide, it is also worth noting that future X-ray missions such as XEUS and 
Generation-X 5 should allow the first investigations of the evolution of the normal early-type galaxy 
X-ray luminosity function with redshift. 

5. Summary 

Using X-ray stacking analyses, we have investigated the X-ray evolution of 539, z « 0.1-0.7 
early- type galaxies located in the E-CDF-S. These galaxies were selected using a combination of 
red-sequence colors and Sersic indices as a part of the COMBO-17 and GEMS surveys (M05). 
We classified our original early-type galaxy sample as the "general sample" and generated an 
additional "faded sample," which was corrected for the passive fading of old stellar populations. 
Using these samples, we analyzed separately optically luminous (Lb ~ 10 10-11 Lb,®) and faint 
(Lb ~ 10 93-10 Lb,®) populations, which are expected to have soft X-ray spectra dominated by 
hot interstellar gas and LMXBs, respectively. Our primary goal was to use stacking analyses to 
measure and constrain the redshift evolution of the average X-ray emission from normal early-type 
galaxies. To achieve this, we used a variety of techniques to identify powerful AGNs, which we 
removed from our stacking analyses. Our key results follows: 

1. We detected 49 early- type galaxies in the X-ray band and classified 32 of these as AGN 
candidates based on their X-ray, optical, and radio properties (see § 2.1.1 for details); the 
remaining 17 X-ray-detected sources had multiwavelength properties consistent with normal 
galaxies. In addition to the 32 X-ray-detected AGN candidates, we identified 13 galaxies 
with AGN-like radio-to-optical flux ratios, which we characterized as potential AGNs. We 
found that the majority of the AGN candidates were coincident with optically luminous early- 



5 For further information regarding the future X-ray missions XEUS and Generation-X, see 
http://www.rssd.esa.int/XEUS/ and http://genx.cfa.harvard.edu, respectively. 
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type hosts. The inferred AGN fraction for our optically luminous galaxies shows evidence 
for evolution with redshift in a manner consistent with the (1 + z) 3 evolution expected from 
other studies of AGN evolution. 

2. When stacking the X-ray counts from our normal optically luminous early- type galaxy sam- 
ples, we found that the X-ray-to-optical mean luminosity ratio, Lx/Lb, stays roughly con- 
stant over the redshift range z ~ 0.0-0.7, which indicates that the X-ray-emitting gas has not 
significantly evolved over the last ~6.3 Gyr (i.e., since z ~ 0.7). Using the data available, 
we found that (Lx/Lb) z =o.7 = [1-0 ± 0.5] and [1.7 ± 0.8] x (Lx/ Lb) z =o for our general and 
faded samples, respectively. We interpret the lack of X-ray evolution of optically luminous 
early-type galaxies to be due to an energy balance between the heating and cooling of the 
hot gas over each cooling time. When assuming that the heating is largely due to transient 
AGN activity, we found that mechanical feedback dominates the heating out to z ~ 0.7 versus 
radiative power, which we inferred to be very poorly coupled to the gas. Furthermore, this 
result suggests that the radiative and mechanical AGN power evolve differently with cosmic 
time. 

3. For our optically faint early- type galaxy samples, we found suggestive evidence that Lx/Lb 
increases with redshift. By z « 0.5, Lx/Lb is measured to be 5.3 ± 4.1 and 5.6 ± 4.1 times 
larger than that measured at z = for our general and faded samples, respectively; however, 
due to poor statistical constraints on the local relation and the undetected AGN population, 
we could not confidently rule out the null hypothesis. We hypothesized that evolution of the 
optically faint early-type galaxy X-ray emission may be due to the evolution of LMXBs in 
galaxies that have recently joined the red-sequence and/or the cooling of hot gas within these 
galaxies. 
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Table 1. X-ray Detected Early- Type Galaxies: Source Properties 
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Table 1 — Continued 
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Note. — Col.(l): Chandra source name. Col. (2): Source redshift as determined by COMBO-17. Col.(3J-(4J: Flux for the 0.5-2.0 keV and 2-8 keV bandpasscs. Col. (5): 
Hardness ratio of the 2—8 kcV and 0.5—2.0 kcV count rates (^2 — 8 fceV /^0.5 — 2.0 kev)- Col. (6): Logarithm of the 0.5—8.0 kcV to .K-band flux ratio. Col. (7): Effective photon 
index Here a value of 1.4 was assumed when photon statistics were too limited to determine accurate values. Col. (8): Logarithm of the rest- frame B-band luminosity. 

Col. (9): Logarithm of the rest-frame 2—8 keV luminosity. Col. (10): Survey field in which each source was identified. For E-CDF-S identifications, the associated field number 
(i.e., 01— 04J indicates the Chandra pointing within which the source was detected (sec Lchmer ct al. 2005b for details). Col. (11): Source notes. Here an "N" denotes normal 
galaxies, an "A" denotes candidate AGNs, an "R" denotes sources with radio detections, and an "S" denotes sources that were detected in the supplementary catalogs (i.e., 
using wavdetect at a false-positive probability threshold of 1 X It) - 5 ). 



Table 2. Stacked Early- Type Normal Galaxies: Basic Properties 
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Note. — Col.(l): Mean rcdshift and standard deviation of the redshift for each stacked sample. Col. (2): Number of sources being stacked. Col. (3): 
Number of stacked sources that were detected individually in the X-ray bandpass. Col. (4): Total vignctting-corrcctcd effective exposure time measured 
from the 0.5-2.0 kcV exposure maps. Col.(5)-(7): Nct 0.5-1.0 kcV, 0.5-2.0 kcV, and 2-8 kcV source counts. Col. (8J-(10) : S/N for the 0.5-1.0 kcV, 
0.5-2.0 kcV, and 2-8 kcV bandpasses. 



Table 3. Stacked Early- Type Normal Galaxies: Mean X-ray Properties 
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0.22 ± 0.09 
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a For the faded samples, quoted S-band luminosities indicate faded (z — 0) luminosities, 1/^,0- 



Note. — Col.(l): Mean rcdshift and standard deviation of the rcdshift for each stacked sample. Col. (2)— (4): Logarithm of the mean 0.5—1.0 keV, 0.5—2.0 keV, and 2—8 keV flux in units 
of ergs cm -2 s — 1 . Col. (5): Logarithm of the 0.5—8.0 keV to K-band flux ratio. Col. (6)— (8): Logarithm of the mean 0.5—1.0 keV, 0.5—2.0 keV, and 2—8 keV rest-frame luminosity. Col. (9): 
Logarithm of the mean rest-frame B-band luminosity. Col. (10): Logarithm of the 0.5—2.0 keV— to— S-band mean luminosity ratio. 
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Fig. 1. — Rest-frame X-ray spectral energy distributions (SEDs) for (a) NGC 1600 (Sivakoff etal. 2004) and 
(b) NGC 4697 (Sarazin etal. 2001). These are normal early-type galaxies with X-ray emission dominated by hot 
interstellar gas (NGC 1600) and discrete LMXBs (NGC 4697) and have X-ray spectra representative of optically 
luminous and faint galaxies, respectively. These SEDs (solid curves) were separated into LMXB (dashed curves) 
and hot gas (dotted curves) components (see § 2.2 for further details). The X-ray emission from NGC 1600 was 
best fit using a power-law (r « 1.8) for the discrete component and a MEKAL plasma (kT w 1 keV; Z w 0.2Zq) 
for the unresolved component. NGC 4697 was best fit by a thermal bremsstrahlung (kT m 5.2 keV) for the LMXB 
component and a MEKAL plasma (kT w 0.3 keV; Z w O.IZq) for the hot interstellar gas. For reference, we have 
shown the energy ranges of our adopted bandpasses redshifted to z = 0.55 and z = 0.39, the median redshifts of our 
optically luminous and faint samples, respectively. We note that for X-ray luminous (optically luminous) early-type 
galaxies within the redshift ranges considered in this study, the 0.5-1.0 keV (SB1) and 0.5-2.0 keV (SB) bandpasses 
primarily trace hot interstellar gas and the 2-8 keV (HB) bandpass is primarily dominated by LMXBs. However, for 
X-ray faint (optically faint) early-type galaxies, all bandpasses are generally dominated by LMXB emission. 
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Fig. 2. — (a) Rest-frame _B-band luminosity versus redshift for our general sample of 539 early-type galaxies. Larger 
circles indicate X-ray-detected AGN candidates (open) and normal galaxies {filled). Sources denoted with crosses are 
(1) in close proximity (<10") to an X-ray-detected source, (2) within the boundaries of extended X-ray sources, (3) at 
large off-axis angles (i.e., >7' from all aimpoints), and/or (4) found to have relatively large radio-to-optical flux ratios 
(see § 2.1.1); these sources have been removed from our stacking analyses (see § 2.2 for details). The shaded bands 
show the luminosity ranges of our optically luminous and faint samples; vertical dotted lines indicate the evenly-spaced 
comoving volume intervals chosen in constructing our stacking samples, (b) Evolved, z = 0, B-band luminosity of 
our sample (see discussion in § 2.1.2), which constitutes our faded sample. Symbols and boundaries are the same 
as in Figure 3a. Using a Chabrier initial mass function (Chabrier 2003), we estimate that the luminosity ranges 
L B ,o = 10 9 ' 3 " 10 and lO 10 " 11 L B ,o correspond roughly to stellar mass ranges of « lO 99 " 106 and « lo 10 6 " 11 - 6 M , 
respectively. 
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Fig. 3. — (a) Soft-band (0.5-2.0 keV) signal-to-noise ratio (S/N) versus inclusion radius (i.e., the maximum off-axis 
angle within which a source is included in stacking) for a fixed aperture size of l'.'5; the solid line indicates the expected 
linear increase in S/N versus inclusion radius, (b) Soft-band S/N versus aperture radius for a fixed inclusion radius 
of 7'.0. The above plots were generated using corresponding optimized stacking parameters (i.e., stacking aperture 
size and inclusion radius), which were determined iteratively (see § 2.2 for details); the optimized parameters are 
indicated by vertical dashed lines in each plot. 
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Fig. 4. — Adaptively smoothed 0.5-2.0 keV image of the combined wl Ms CDF-S and ^250 ks E-CDF-S. Positions 
of stacked sources from our faded sample are shown as open circles (ml Ms CDF-S) and open diamonds (~250 ks 
E-CDF-S); filled circles represent galaxy positions that have been stacked using both the wl Ms CDF-S and «250 ks 
E-CDF-S observations (see § 2.2 for additional details). Aim points of each Chandra observation are indicated as 
plus signs, and the surrounding 7^0 inclusion radii are indicated with dashed circles. The apparent lack of sources 
in the north-eastern corner (i.e., the upper left-hand corner) of the image is partially due to missing HST coverage 
from the GEMS imaging. For reference, we have outlined the «160 arcmin 2 GOODS-S region (solid rotated rectangle; 
Giavalisco etal. 2004). 
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Fig. 5. — Optical 72-band magnitude versus 0.5-8.0 keV flux for X-ray-detected early-type galaxies in our sample. 
AGN candidates and normal galaxies are plotted as open symbols (circles and stars) and filled circles, respectively; 
sources with upper limits were detected in either the 0.5-2.0 keV or 2-8 keV bandpasses. Diagonal dotted lines 
represent lines of constant X-ray-to-optical flux ratio (i.e., log/0.5-8.0 kcv/fn =+1, — 1, and —2); luminous AGNs 
generally have log fo. 5-8.0 kov / fa £ — 1 (shaded region). AGN-candidates were classified following the three criteria 
discussed in § 2.1.1. The sources marked with five-pointed stars are radio-detected AGN candidates, including the 
FR II source CXOECDFS J033228. 81-274355. 6. Furthermore, mean i?-band magnitudes and 0.5-8.0 keV fluxes for 
our stacked optically luminous (large filled squares) and faint (large filled triangles) faded samples are plotted with 
varying grayscale levels to indicate redshift; darker levels indicate larger redshift values. 
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Fig. 6. — Cumulative X-ray detected AGN fraction, fc, as a function of 2-8 keV luminosity for our optically 
luminous faded samples with z m 0.10-0.55 (filled circles) and z w 0.55-0.70 (open diamonds). Each data point with 
printed fractions represents the number of AGN candidates detected with 2-8 keV luminosity of L2-8 kcv or greater 
divided by the number of early-type galaxies with X-ray coverage sufficient to detect an AGN of L2-8 k c v (see § 3.1 
for further details). 
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Fig. 7. — Stacked, adaptively-smoothed 0.5-2.0 keV images of our early-type, faded samples for each redshift bin. 
The top panels show our optically luminous {Lb,o ~ 10 10 ~ n Lb,®) samples, and the bottom panels show our optically 
faint (Lb,o ~ 10 9 3-10 Lb,q) samples. These images were generated using the CIA0 tool CSM00TH with a minimum 
significance of 2.5a. The images are «15" (~30.5 pixels) per side, and each pixel is 0'.'492. Our circular stacking 
aperture of radius l'.'5 is shown in each image centered on the optical centroid of our stacked sources. Additional 
sample information, including the number of stacked galaxies, total Chandra exposure, signal-to-noise ratio (S/N), 
and logarithm of the mean z — B-band luminosity {Lb,o) are annotated on each smoothed image. 
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Fig. 8. — Logarithm of the 0.5-2.0 keV luminosity (Lx) versus the logarithm of the £?-band luminosity, Lb, for 
D < 70 kpc local ellipticals with the average properties of our faded samples plotted. Small symbols and upper limits 
are sources from the OS01 local sample, and the different symbols correspond to field galaxies (circles), central-cluster 
galaxies (upward-pointing open triangles), brightest-group galaxies (downward-pointing open triangles), and AGNs 
(crosses). The best-fit relations for the luminosity intervals Lb ~ 10 1() Lb,q and Lb ~ 10 10 Lb,q are shown as solid 
and dashed lines, respectively; the shaded region shows the expected discrete-source contribution (Kim & Fabbiano 
2004; sec also the discussion in § 3.2.1). The local ellipticals M32, M87, NGC 1399, NGC 1600, NGC 4697, and 
NGC 5102 have been marked, for reference. Large shaded symbols represent our optically luminous (squares) and 
faint (triangles) faded samples and have grayscale levels corresponding to the mean redshift of each sample, such 
that darker shading represents higher redshifts. Error bars in Lx represent la errors on the mean. We quote £?-band 
luminosities as Lb,o to illustrate any potential evolution of the mean X-ray luminosities. 
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Fig. 9. — Logarithm of the 0.5-2.0 keV luminosity Lx versus redshift for our general (a; see Figure 2a) and 
faded (b; see Figure 2b) samples; filled and open circles indicate X-ray-detected normal galaxies and AGN candi- 
dates, respectively. Black squares and triangles with error bars show the stacking results for our optically luminous 
(L B ~ 10 10 ~ n ergs s _1 ) and faint (L B ~ 10 9,3 ~ 10 ergs s _1 ) samples, respectively. Error bars in redshift represent 
the standard deviation of the redshift for sources in each stacked sample. For comparison, we have plotted the corre- 
sponding mean X-ray luminosities (and errors on the means computed with ASURV) of normal early-type galaxies from 
the OS01 local sample (gray filled square and triangle). All mean values (ours and those of OS01) were calculated 
after excluding AGNs, central-cluster galaxies, and brightest-group galaxies, and should reflect the average properties 
of isolated field early-type galaxies. The solid curve illustrates the median X-ray detection limit (using the median 
sensitivity limit of x 10 -16 ergs cm~ 2 s _1 for our total sample). For reference, the X-ray luminosities of the 

local ellipticals M32, M87, NGC 1399, NGC 1600, NGC 4697, and NGC 5102 have been indicated. 
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Fig. 10. — Logarithm of the X-ray-to-B-band mean luminosity ratio (Lx/Ls) versus redshift for our general (a; 
see Figure 2a) and faded (b; see Figure 2b) samples. Symbols have the same meaning as in Figure 9; the dotted line 
and shaded region represent the expected local discrete-source contribution and its dispersion (from Kim & Fabbiano 
2004). For the faded sample, we used Lb,o when computing the X-ray-to-B-band mean luminosity ratios. For 
reference, the local ellipticals M32, M87, NGC 1399, NGC 1600, NGC 4697, and NGC 5102 have been plotted. 
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Fig. 11. — Residuals to the local best-fit relation, log Lx — a log Lb — j3 for our general (a; see Figure 2a) and faded 
(b; see Figure 2b) samples. For optically luminous galaxies (Lb ~ 10 10 " 11 Lb,®) we used a = 2.61 and fi = 12.77, 
and for optically faint galaxies (Lb ~ 10 9,3 ~ 1() Lb,q), we adopted a = 1.05 and /? = 29.36. Symbols have the same 
meaning as in Figure 9, and the redshifts of the local samples have been offset from z = for viewing ease. The 
dotted horizontal line indicates the zero residual. For reference, the local ellipticals M32, M87, NGC 1399, NGC 1600, 
NGC 4697, and NGC 5102 have been plotted. 
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Fig. 12. — (a) Cumulative AGN fraction (i.e., the fraction of galaxies harboring an AGN with a 2-8 keV luminosity 
of 1/2-8 keV or greater), fc, versus logZ/2-8 k c v for our optically luminous faded samples. The cumulative AGN 
fraction computed over the entire redshift range, {fc}z, is indicated as filled circles with la error bars, and the 
thick solid curve represents our best-fit quadratic relation to the data (see § 3.2.2 for details). This relation was 
used to estimate fc for our optically luminous faded samples (annotated curves). For reference, we have plotted 
the observed AGN fraction for 2 « 2.5 DRGs (Rubin et al. 2004; open diamond) and our model extrapolated out to 
2 = 2.5 (dot-dashed curve), (b) Differential AGN fractions (i.e., the fraction of galaxies harboring an AGN in discrete 
bins of width AlogL2-8 k c v = 0.5), fo, versus logL2-8 keV- (c) Fraction of early-type galaxies for which we could 
not have detected an AGN with a 2-8 keV luminosity of L2-8 kov, /s, versus logZ/2-8 keV- (d) Fraction of galaxies 
harboring AGNs in our optically luminous faded samples that would remain undetected due to sensitivity limitations, 
fu = Jd x /s, versus logZ/2-8 k c v; these galaxies would not have been removed from our stacking analyses. 
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Fig. 13. — Mean 2-8 keV AGN luminosity per early-type galaxy, L2-8 k c v,AGN, versus redshift for our opti- 
cally luminous faded samples (filled squares); these luminosities were derived following equation 8. Error bars on 
L2-8 kcV.AGN are la errors, which were derived by propogating errors on the fit to {fc)z (solid curve in Fig. 12) 
through to equation 8. For comparison, we have plotted the Brand et al. (2005) stacking constraints for their early- 
type galaxy samples (open diamonds), which have mean i?-band absolute magnitudes similar to those of our optically 
luminous faded samples. The dotted curve represents the (1 + z) 3 model fit to our data. 



